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The transfer matrix DMRG method for one dimensional quantum lattice systems has been devel- 
oped by considering the symmetry property of the transfer matrix and introducing the asymmetric 
reduced density matrix. We have evaluated a number of thermodynamic quantities of the anisotropic 
spin-1/2 Heisenberg model using this method and found that the results agree very accurately with 
the exact ones. The relative errors for the spin susceptibility are less than 10"^ down to T = O.OIJ 
with 80 states kept. 
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The density matrix renormahzation group (DMRG) j|] 
is a powerful numerical method for studying the ground 
and low-lying state properties of low dimensional lattice 
models. It has been applied successfully to a number of 
strongly correlated systems at zero temperature in one di- 
mension (Id) (2|H] as well as two dimensions (2d) 
In 1995, Nishino pointed out that as the partition func- 
tion of a 2d classical system is determined only by the 
maximum eigenvalue of a transfer matrix in the thermo- 
dynamic limit, the DMRG idea can be extended to find 
this extreme eigenvalue and the corresponding eigenstate. 
Subsequently one can study thermodynamic properties. 
He calculated the specific heat of the 2d Ising model us- 
ing this so called transfer matrix DMRG method and 
found that the result agrees very accurately with the ex- 
act solution jio). Recently, Bursill, Xiang and Gehring 
have developed a transfer matrix DMRG algorithm for 
Id quantum systems They tested the method on 

the dimerized spin 1/2 XY model and obtained encour- 
aging results. However, a proper implementation of the 
transfer matrix DMRG for studying the thermodynamic 
properties of Id quantum systems, especially at low tem- 
perature, remains challenging. 

In this paper, we report our recent progress on the 
development of the transfer matrix DMRG method for 
Id quantum systems. We have considered the symmetry 
properties of the quantum transfer matrix [T^JT^ and in- 
troduced an asymmetric reduced density matrix which 
optimizes truncated basis states. This study allows us to 
achieve significantly accurate results at low temperature 
and greatly enhances the applicability of the method. 
For the S = 1/2 Heisenberg antiferromagnetic chain, 
we found that the relative error for the spin suscepti- 
bility is of the order of 10~^ down to the temperature 
T — O.OIJ (J is the exchange constant) with 80 states 
kept. This size of error and the value of the temperature 
are smaller than those for typical Quantum Monte Carlo 
results as well as the thermodynamic DMRG results . 
The thermodynamic DMRG method |l^ cannot treat ac- 
curately a physical system where the correlated length 



diverges (for example the low temperature region of the 
S=l/2 Heisenberg model) because of the finte size effect. 
However, there is no such problem in the transfer matrix 
DMRG method since this method deals directly with an 
infinite lattice system. 

We demonstrate our method using the Id anisotropic 
spin-1/2 Heisenberg antiferromagnetic model: 
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We shall set J = 1. To apply the DMRG idea, we use 
the Trotter formula to decompose the partition function. 
H is separated into two parts, Ho and He, containing 
those terms with i being odd or even, respectively. The 
partition function is represented in terms of the quantum 
transfer matrix Tm [l2lE3|: 
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where e = (3/M, (3 = 1/T and M is the Trotter number. 
The elements of the asymmetric matrix Tm are deter- 
mined by the product of 2M local transfer matrices 
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with T{crlMcrlMj^^\alj^,ja\M+i) - '^(^'Ia/'^i l"2Af 
imposing periodic boundary conditions in the Trot- 
ter direction. The local transfer matrix is given by 

^ =(.;+^,4+i|exp(-e/.,)l4,^r^') 
where cr\ = {—'^y^'^s]. and \s\) is an eigenstate of S'f: 

^iWk)—^kWk)- The basis |cr^) ® \'^k+i) is used to repre- 
sent r and to construct the corresponding Trotter space. 
The superscripts and subscripts in Tm and t represent 
the coordinates of spins in the real and Trotter space, 
respectively. 

The local Hamiltonian hi conserves the total spin 
at sites i and i + l, i.e. s\ + s^"^^ = s^.^^ -I- s^^\. 
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In the Trotter space this conservation law can be ex- 



pressed as al + a 
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crj.+^ + al:^j^. This means that 



T(a^"'"^cr^'^\|CT^cr^_j,j) is block diagonal for each value of 
(7^ + o'fc+i- It turns out that the total sum of a^. at site 

i is conserved in Tm, i.e. J^k'^l ^ Efe '^fe^^ "^m 
is block diagonal according to the value of '^j,cr\- For 
Eq. (1), it can be further proved that the maximum 
eigenvalue of Tm occurs in the =0 subblock [p^ . 

Therefore only the 'Ylk'^k ~ ^ subblock in Tm is consid- 
ered in our DMRG iterations. This consideration allows 
to keep more basis states in the truncation of the basis 
set and to save computer CPU time. 

In the limit iV — > oo, one can study the thermody- 
namic properties using the maximum eigenvalue A and 
corresponding left ("0^1 and right \^^) eigenvectors of 
the transfer matrix T/v/. The free energy is determined 
purely by A, F = — ^InA. From the derivatives of F 
one can in principle calculate the internal energy U , the 
magnetization M^, the specific heat Cy, the spin suscep- 
tibility X) other quantities. However, as it is difficult 
to evaluate accurately a derivative of a function in numer- 
ical calculations, we find that it is better to evaluate U 
and Mz directly from and and then calculate 

Ct} and X from the derivative of U and Mz , respectively. 
For instance, U^{H)T/N^{hi)T^{iP^\Th^\i>^)/\, where 
( )t is the thermal average with respect to the thermody- 



namic density matrix pth=exp{— PH) / Z and {i/j^lip^)— 1 
is assumed hereafter. The definition of 7^^ is similar to 
that of Tm subject to the decomposition. Its matrix el- 
ements can be obtained from the right hand side of Eq. 
(|) by replacing T{alal\a\a\) with Th^{(jlal\a\a\) = 
(sjSjIft-i exp(— t/ii)|s2Si). Similarly, one can find out 
the relation between the magnetization Mz=(^- S^t/N 
and the maximum eigenvectors of Tm- 

In our calculation, we fix e and increase the chain 
length 2M. For each M, the temperature T — 1/eM. 
As M is small, one can find A, {tp^l and 1-0^} exactly. 
For large M we extend the DMRG idea to approximately 
but accurately find A, {ip^l and for a periodic time- 
slice chain. Figure 1 shows two configurations of the su- 
perblocks of Tm- The superblock consists of two blocks, 
which we call renormalized blocks, in the dashed frames 
and two time slices. The system contains a renormalized 
block and one slice. The rest is thus its environment. We 
use Us and Ue to label the basis states of the renormalized 
blocks in the system and the environment, respectively. 
The states of two time slices are represented by cti and 
CT2. The elements of the right transfer matrix is denoted 
by To{a'{,n'^,a2]ai,ns,a2) or T; ((7^, n^, cr^'; 0"^', n^, (72) if 
M = odd, or even. The left transfer matrix can be ob- 
tained by transposing the right one. Therefore the su- 
perblocks Tm is given by: 
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To grow the chain, we have the following recursive relations: 

T (ct'i ,n'^,a2;a'{,ns,a2) ^^t{(t[,<t'\ a'{ , a" ) % (a" , n(, , g'{ ; a, , 02 ) 
To {a'i ,n'^,a2;ai,hs,(J2) = y^T{a'{,(T"\(Ti,a)T{a',n'^, (J2 ; cr" , , 0-2 ) 
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which enlarge the renormalized blocks by one slice with \hs) = \<t) ® \ns) as the corresponding Trotter space and \a) 
as the states of a spin added. Initially, for 2M — 4, X=(cri, cr', CTj ; cr", cr, 0-2) = X^o-" '''{^'n'^'Wi^ ^'^''iWi '^2)- As 

the number of states in \hs) exceeds m, Te{a[, n'g,a2; a i,hs, (72) and 7^(cr", rig, cr2 ; cti, n^, CT2) are renormalized by: 



T{a[,n'^, cr^,'; a'{, n^, (72) = ^ 0\n'^,n'JT{a[,h'^, a!^; cr'l, n^, 0-2)0'' (n^, n^), 
To{a'{, n'^, cr^'; cti, n^, era) = ^ 0\n'^,n'^)To{a'{, n'^^cr'^; cti, ri^, CT2)0''(ns, n^), 
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where the transformation matrices O'''' are constructed by m truncated basis states from a reduced density matrix 
discussed below. Other operators such as hi can be renormalized by Eq. (6) with th^ instead of r in Eq. (5). 



We compute the maximum eigenvalue. A, and the cor- 
responding right eigenvector, \tp^), of Tm using a pro- 
jection method. Iterating \iPk) = Tm\iJjk-i), we reach 
= \iPk) and Tm\iI^^) = X\ip-'^') for sufficient large 
K- j-^o) is an arbitrary trial vector which is not or- 



thogonal to In our calculations, we find that the 

value of K needed for producing an eigenvalue with a 
relative error less than 10~^^ is generally less than 20, 
but it increases with increasing M. The left eigenvec- 
tor \^p^) can be calculated similarly. However, for sys- 
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terns as we study here, the wave function of (■0^1 can 
be directly read out from the wave function of \%p^): 
■>p^{ns,cr2,cri,ne) = ip^{ne,a2,o-i,ns) by constructing 
the superblocks with a reflection symmetry as involved 
in Eq. (4). 

A density matrix for the whole system (i.e. superblock) 
can be defined as p = T^^^ /TtT^^^ . This is a general- 
ization of the thermodynamic density matrix pth in the 
Trotter space. We form the reduced density matrix for 
the augmented renormalized block by performing a par- 
tial trace on p for the states of the environment 
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In the thermodynamic limit, ps — Tr„^o-2 thus 
the matrix element of ps is given by 



(8) 



with \n) = \a) ®\n). ps is an asymmetric matrix since 
{tp^\ ^ dV'^})^- The eigenvalue of ps gives the prob- 
ability of the corresponding eigenstates onto which the 
system is projected as the response to its environment. 

which is used to define the density matrix for 
the augmented system block in Ref. jlll is not a true pro- 
jection operator for the maximum eigenvectors of Tjf.) 
The transformation matrices O''"" in Eq. (6) are thus 
built up by using m left or right eigenvectors of ps cor- 
responding to m most probable eigenvalues. 

Systematic errors come from two sources. One is the 
finiteness of e, and the other is the truncation of basis 
set in the DMRG iterations. The first type of error is 
generally very small and in principle it can be further 
reduced by doing an extrapolation with respect to [?] 
The error due to the truncation is difficult to estimate. 
A lower bound for this type of error is given by the trun- 
cation error, p„i = 1 ^ YIT '^i^ where Wi {i = 1,- • •, m) 
are the m largest eigenvalues of ps. We found that Pm 
is generally less than 10~^ when m = 16 and decreases 
rapidly with increasing m for the spin 1/2 system. 

Figure 2 shows the results for the specific heat Cv{T) 
down to T = 0.02 with to = 80 and e = 0.05 for A = 0, 1 
cases. Cv is obtained from the first derivative of U. For 
the XY model (A=0), we find that the relative errors 
are less than 10~^ down to T = 0.02 for U and less than 
lO^'^ down to T = 0.03 for C„ compared with the exact 
results. For the isotropic antiferromagnetic Hciscnberg 
model (A=l), the precision of the results is similar. The 
maximum value of is 0.3515 at T = 0.47. At low tem- 
perature Cy varies linearly with T. The coefficient of the 
T term is shown to be 26'/(3sin6') with 9 = coa'^ A Q. 
By fitting our results with a polynomial up to 7th order 
in T for 0.03 ~ T < 0.1, we found that the coefficients 
of the linear terms are 1.041 and 0.665 for A = and 1, 



respectively. The difference between our results and the 
exact ones for the linear coefficients is less than 1%. 

For comparison, we also calculated U and Cy with a 
symmetric density matrix as defined in Ref. [pl| . At high 
temperature, the results obtained with a symmetric den- 
sity matrix agree well with those obtained with an asym- 
metric density matrix. However, at low temperature we 
found that the results obtained with an asymmetric den- 
sity matrix are more accurate than those obtained with a 
symmetric density matrix (Figure 2). The relative errors 
for U and Cy obtained with the symmetric density matrix 
are generally larger than 10^^ at low temperature. 

Figure 3 shows our results for the spin susceptibility 
X(r) down to T = 0.01 with to = 80 and e = 0.05 for 
A = 0, 1 cases. x{T) is obtained from the first order 
derivative of M^, which is equal to M^iT, B)/B for suf- 
ficient small magnetic field since Mz{T,0) — 0. In our 
calculations, B = 0.003 is used. For both cases, the rel- 
ative error is less than 10~^ down to T = 0.01. We note 
that the results of spin susceptibility are generally more 
accurate than those of the specific heat. In the inset, we 
compare numerical results for A = 1 with m = 32 and 
80 to the exact results in the low temperature regime. 
For TO = 32 the relative error is of the order of lO"'^ 
at T = 0.01. Our results are systematically better than 
those obtained by Moukouri and Caron with the thermo- 
dynamic DMRG method ||l|l. 

Our computations were performed on DEC Alpha sta- 
tions. It takes about 14000 seconds on a 175MIIz station 
to generate a superblock size of 2M=4000 for m — 32. 

In conclusion, the quantum transfer matrix DMRG 
method with asymmetric density matrices is developed. 
We have calculated a number of thermodynamic quan- 
tities for the anisotropic Heisenberg antiferromagnetic 
spin-1/2 model and found the results agree very accu- 
rately with exact ones. Our investigation shows that the 
transfer matrix DMRG is a very promising method for 
studying thermodynamic properties of Id quantum sys- 
tems. 
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FIG. 1. Configurations of the superblocks: (a) M = odd 
and (b) M = even. The left and right transfer matrices are 
connected, via the summation over states tr" and (72, to form 
a periodic time-slice chain for Tm- 

FIG. 2. The specific heat C„(T) for both A = and 
A = 1. The solid curve is the exact results. Circles and 
pluses are the transfer matrix DMRG results with the asym- 
metric and symmetric density matrix, respectively, e = 0.05 
and m = 80 are used in the transfer matrix DMRG calcula- 
tions. Inset: a polynomial fit for the low temperature C„. 

FIG. 3. The spin susceptibility x{T) for both A = and 

A = 1. The solid curves and circles, respectively, are the exact 
results and the transfer matrix DMRG results with e = 0.05 
and m = 80. Inset compares the transfer matrix DMRG re- 
sults for m = 32 (diamonds) and 80 (circles) with the Bethe 
ansatz results in the low temperature regime for A = 1. 
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Wang et al Fig. 1(a) 



Wang et al Fig. 1(b) 
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